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1.  INTRODUCTION 


This  report  describes  work  completed  and  in  progress  on 
the  development  of  a  fortran  computer  code  FASCOIi?  to  calculate 
atmospheric  transmittance  and/or  radiance  for  any  given  path  at 

high  spectral  resolution.  The  present  code  is  based  on  FASCODE^ 

2  '-i 

and  FASC0D1  computer  codes  that  use  the  AFGL  line  compilations' 

and  an  approximation  to  the  Voigt  line  profile  based  on  the 

4 

HIRACC  algorithm.  All  aspects  of  this  code  that  are  indepen¬ 
dent  of  spectral  resolution  are  to  be  the  same  as  those  found 
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in  the  moderate  spectral  resolution  code  LOWTRANU  ,  once  it 
has  been  modified  to  accept  the  improved  geometry  routine  dis¬ 
cussed  here.  Most  importantly,  the  two  codes  have  a  common  set 
of  model  (gas  and  aerosol)  atmospheres  and  are  run  by  similar 
sets  of  control  cards. 

Section  2  describes  how  the  code  has  been  modified  to  cal¬ 
culate  atmospheric  paths  for  a  spherical,  refracting  medium 
using  a  new  AFGL  geometry  routine^  along  with  the  six  represen¬ 
tative  model  atmospheres  found  in  the  LOWTRAN  codes.  Following 
this,  the  code  has  been  changed  to  include  zeroth-order  scat¬ 
tering  and  absorption  due  to  the  aerosol  models  used  in  the 

O 

LOWTRAN  5  code  .  These  last  changes  are  detailed  in  Section  3 
and  the  resulting  interim  code  I  documented  in  the  Appendix. 

FASCODE's  radiative  transfer  assumes  the  atmospheric  layers 
are  in  local  thermodynamic  equilibrium  (LTF.).  For  higher  alti¬ 
tudes,  the  photon-driven  chemistry  and  reduced  molecular  colli- 
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sion  rate  cause  this  assumption  to  be  invalid.  Degges  has 
modelled  the  kinetics  of  this  phenomena  and  produced  selected 
single-transition  population  densities  for  given  sun  angles 
(season  and  time  of  day).  Section  4  discusses  how  the  results 
from;  Degges  computations  can  be  used  to  include  upper'  atmospher¬ 
ic  non-LTE  effects  In  FASCODE's  transmittance  and  radiance  cal¬ 
culations  . 

The  Voigt  line  shape  model  predicts  a  Lorenzian  profile 


for  the  line  wings.  It  has  been  established  that  this  is  in¬ 
correct  for  CO  .  The  far  wings  are  sub-Lorentzian .  This  error 
becomes  important  near  band  heads  located  next  to  window  regions 
in  the  atmosphere  where  the  background  continuum  constitutes  the 
only  absorption.  The  focus  in  the  present  work  is  the  CO^  con¬ 
tinuum.  Section  5  outlines  this  problem  and  discusses  possible 
ways  the  line  profiles  can  be  modified  in  the  far  wing  to  ac¬ 
count  for  this  aspect  of  atmospheric  radiative  transfer. 

The  inhomogeneous  model  atmospheres  are  approximated  in 
FA3C0DE  by  a  user-specified  number  of  homogeneous  layers  to  sim¬ 
plify  the  radiative  transfer.  Consequently,  line  strengths  and 
half  widths  are  evaluated  at  a  defined  average  therodynamic 
state  in  each  layer.  This  approach  produces  errors  that  can 
only  be  decreased  by  increasing  the  number  of  layers  and  there¬ 
fore  increasing  the  computational  time.  The  proposed  study  out¬ 
lined  in  Section  6  is  to  develop  some  knowledge  concerning  how 
layering  produces  errors  so  that  they  can  be  controlled.  The 
goal  of  this  work  is  to  produce  guidelines  for  selecting  layer 
boundaries  to  minimize  the  error  while  keeping  the  number-  of 
layers  and  thus  computational  time  acceptable. 
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2.  OPTICAL  PATHS  AND  ATMOSPHERIC  MODELS 

Program  ATMPTH^  and  its  subroutines,  as  supplied  by  AEGL, 
eliminate  the  need  for  offline  programs  such  as  LAYER**  and 
DRIVER  in  the  preparation  of  the  atmospheric  information  for 
PA3C0D1 .  The  program  is  capable  of  determining,  optical  paths 
through  the  atmosphere,  and  includes  the  effects  due,-  to  spheri¬ 
cal  geometry  and  refraction.  ATMPTII  calculates  the  molecular 
density-weighted  average  temperature  and  pressure,  along  with 
the  Integrated  molecular  absorber  amounts  in  each  of  the  homo¬ 
geneous  layers  used  to  approximate  the  atmosphere.  The  user  can 
choose  between  six  representative  atmospheres,  or  supply  his 
own.  The  number  and  placement  of  layer  boundaries  is  also  spec¬ 
ified  by  the  user.  The  goal  is  to  incorporate  ATMPTH  as  a  sub¬ 
routine  to  FASC0D1,  eliminating  the  need  for  off-line  programs. 
Description  of  the  ATMPTH  Program;  A  brief  description  of  the 
ATMPTH  program  will  be  presented  here  for  those  users  who  may 
require  a  more  complete  understanding  of  the  overall  calculation 
procedure.  The  first  control  card  specifies  the  choice  of  model 
atmosphere.  There  are  six  built-in  representative  atmospheres 
which  are  specified  by  the  pressure,  temperature ,  and  molecular 
densities  of  H20,  CO^,  0^,  N^O,  CO,  CH^,  0^ ,  and  at  3A  al¬ 
titudes.  These  model  atmospheres  are  essentially  the  same  as 
those  in  LOWTRAN  ,  except  that  the  densities  are  in  different 
units.  One  of  the  six  model  atmospheres  stored  in  subroutine 
MDLATM,  or  a  non-standard  user-supplied  atmosphere  as  read  by 
subroutine  NSMDL,  is  selected  and  stored  in  the  common  block 
/MDATA/. 

The  second  control  card  must  contain  three  of  the  following 
five  pash  parameters:  HI,  H2,  ANGLE,  RANGE,  and  BETA.  See  the 
LOWTRAN  report  for  allowable  combinations  of  these  parameters. 
Next;  the  desired  PASCODE  layer  boundaries  are  read  in.  These 
need  not  correspond  to  the  altitudes  contained  in  the  model  at¬ 
mospheres  or  include  HI  or  H2,  the  altitudes  of  the  path  end- 
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points.  The  only  restriction  on  the  choice  of  boundaries  is 
that  they  be  chosen  such  that  the  ratio  of  the  mean  half-widths 
in  adjacent  layers  is  not  greater  than  2  to  l\  A  new  atmos¬ 
pheric  profile  containing  both  the  previous  altitudes  and  the 
desired  FASCODE  boundaries  is  formed.  Pressures,  temperatures, 
and  molecular  densities  at  the  FASCODE  boundaries  are  found  by 
interpolation  of  the  model  atmosphere. 

Subroutine  UMTHY  is  called  to  calculate  the  refracted  path 
through  the  atmosphere,  and  integrate  the  amounts  of  the  absorb¬ 
er's  along  that  path.  Subroutine  RFRPTH  (which  drives  the  inte¬ 
gration)  requires  as  input  the  set  111,  ANGLE,  H2  anu  LEN.  These 
are  determined  by  CMTRY  using  the  three  path  parameters  supplied. 
A  final  atmospheric  profile,  generated  by  RERPTH,  starts  at  HMIN 
(tangent  height  if  any)  and  goes  up  the  MAX  (HI,  H2)  (the  larger 
of  111  and  112).  The  profile  is  specified  at  HI,  H2  and  all  of 
the  model  atmosphere  altitudes  and  FASCODE  layer  boundaries  be¬ 
tween  1IMJN  and  MAX  (ill,  H2). 

The  determination  .of  the  refracted  path,  integration  of 
the  absorber  amounts,  and  formation  of  the  homogeneous  layers 
is  done  using  as  layer  boundaries  all  of  the  altitudes  in  the 
final  profile.  Once  the  path  and  integrated  amounts  are  known, 
layers  arc  merged  in  ATMPTH  to  form  the  desired  FASCODE  layers 
defined  by  the  boundaries  HMIN,  HI,  H2,  and  the  FASCODE  bound¬ 
aries  between  HMIN  and  MAX  (Hi,  H2).  Finally,  the  layers  are 
compared.  If  the  ratio  of  the  geometric  thicknesses  of  two 
adjacent  layers  is  greater  -than  10  or  less  than  .1,  the  two 
layers  are  combined  to  form  a  single  layer  with  a  properly 
weighted  average  pressure  and  temperature  and  combined  absorber 
amo  unis . 

Note  that  the  calculation  of  the  refracted  path  and  the 
absorber  amounts  is  done  using,  the  fivodel  atmosphere  altitudes 
as  additional  layer  boundaries.  Only  after  the  calculation  pro¬ 
cedure  is  complete  is  the  number  of  layers  reduced.  This  is 
done  so  that  model  atmosphere  data  is  not  degraded  in  the  process 
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of  performing  the  path  integral:?.  A  general  flow  diagram  of 
ATM PTH  is  given  in  Figure  1. 

ATMPTH  As  A  Subroutine  Of  FASCODI :  ATM  1  Til  and  its  subroutines, 
nave  been  added  to  FAS00D1  as  subrout.  J  lies .  ATM  IT!  I  is  called 
only  once  to  generate  the  layer  pressures  and  tempera  turns ,  in¬ 
tegrated  amounts  of  molecules  and  aerosol:;,  and  the  mean  half- 
widths  associated  with  each  of  the  FA.Ua )I)K  layers.  Tills  infor¬ 
mation  is  stored  in  arrays  in  common  block  /OUTPUT/,  and  is 
passed  to  the  FA3C0DE  subroutine  PATH.  PATH  is  called  as  before 
in  FASCODE,  once  for  each  layer1,  except  that  instead  of  reading 
the  layer  data  from  a  file  generated  by  an  offline  program,  the 
data  is  assigned  using  the  array  values  in  /OUTPUT/  and  /PATH/  a 

PAVE  =  PBAR(L) 

TAVE  =  TBAR(L) 

etc  . 

Due  to  the  number  of  times  that  they  are  referenced  and 
central  memory  considerations,  simple  variables  are  preferred 
over  subscripted  variables  for  use  in  the  calculation  procedures 

The  following  gives  the  sequence  and  format  of  the  ATMPTH 
control  cards.  These  now  appear  as  a  portion  of  the  FASC0D2 
control  cards  and  will  be  mentioneu  again,  but  only  briefly,  in 
the  FA3C0D2  users  guide. 

Three  cards  control  the  operation  of  the  program  while 
other  cards  must  be  read  in  to  define  the  FA PC ODE  boundary  lev¬ 
els,  and  possibly  to  define  non-standard  conditions.  The  for¬ 
mats  of  these  cards  and  the  meaning  of  the  parameters  are  de¬ 
scribed  as  follows: 

CARD  1:  MODEL,  TTY PE,  I IN,  IMOD,  KM AX,  RE 
(515,  5X,  F10.H) 

MODEL  =  0:  USER  SUPPLIED  HORIZONTAL  PATH  PARAMETERS 

1:  TROPICAL  MODEL  ATMOSPHERE 

2:  MIDLATITUDE  SUMMER 


V 

1  1 1 
\ 


L.  CALL  BY  FAS  COD l 

2.  MODEL  ATMOSPHERE  NUMBER 

3.  ATMOSPHERIC  PROFILE 

4.  IMOD,  CALL  FOR  NON-STANDARD  MODEL 
3.  PATH  PARAMETERS 

6.  HI, H2, ANGLE,  LEN 

7.  SINGLE  LAYER  DATA 

8.  INTEGRATED  AMOUNTS  BY  LAYER 

9.  INTEGRATED  AMOUNTS 
10-  BETA,  RANGE,  PHI 

II.  INTEGRATED  AMOUNTS  FOR  FASCOD1  LAYERS 


GMTRY 


3,6 

N 

k  9,10 

RFRPTH 

7 

K  8 

\ 

LAYER 

i 


FIGURE  1  ATMPTH  FLOW  DIAGRAM 


3:  Ml  DLATITUDK  WINTER 

4 :  SUBARCTIC  SUMMER 

5:  SUBARCTIC  WINTER 

6:  U.S.  STANDARD,  1962 

7:  USER  SUPPLIED  ATMOSPHERIC  PROFILE 


ITYPE 


1:  HORIZONTAL  PATH  (CONSTANT  PRESSURE) 

2:  SLANT  PATH  FROM  HI  TO  HP 

3:  SLANT  PATH  FROM  111  TO  SPACE  (100  KM) 


I IN  NUMBER  OF  BOUNDARY  ALTITUDES  FOR  THE  FA3C0C1 
LAYERS  (REQUIRED  FOR  ITYPE  «  2  OR  3) 

IMOD  NUMBER  OF  BOUNDARY  ALTITUDES  FOR  USER 
SUPPLIED  ATMOSPHERIC  PROFILE  (MODEL.  7) 
DEFAULT  =  34 


KM AX  NUMBER  OF  MOLECULAR  SPECIES  FOR  WHICH 

THE  AMOUNTS  aRE  TO  BE  CALCULATED.  DEFAULT  =  8 


RE  RADIUS  OF  THE  EARTH.  DEFAULTS: 

MODEL  =1,  RE  =  6378.39  KM 

2, 3,6,7,  RE  =  6371.23  KM 
4,5,  RE  =  6356.91  KM 

The  formats  for  control  cards  2  and  3  are  different  de¬ 
pending  on  whether  the  path  is  horizontal  (ITYPE  =  1)  or  slant 
(ITYPE  =-  2  or  3) 


For  a  slant  path: 

CARD  2:  HI,  H2 ,  ANGLE,  RANGE,  BETA,  LEN 
(5F10,  4,  15) 

HI  ALTITUDE  OF  THE  OBSERVER  OR  RECEIVER  (KM) 

H2  ALTITUDE  OF  THE  OTHER  ENDPOINT  OF  THE  PATH  (KM) 

ANGLE  ZENITH  ANGLE  AT  HI  (DEGREES) 

RANGE  LENGTH  OF  THE  PATH  FROM  HI  to  112  (KM) 

BETA  EARTH  CENTERED  ANGLE  FOR  THE  PATH  HI  TO  H2  (DEG) 
LEN  =  0,  SHORT  PATH;  =  1,  LONG  PATH  THROUGH  A  TANGENT 
HEIGHT.  LEN  IS  USED  ONLY  WHEN  ANGLE  IS  GT  90.0 
AND  HI  IS  GT  H2 .  DEFAULT  =  0. 


Only  three  of  the  first 
example.  Hi,  112,  ANGLE, 
See  the  comments  in  the 
for  more  details  on  the 


five  parameters  need  be  specified;  for 
or  HI,  H2;  BETA,  or  HI,  ANGLE,  RANGE, 
subroutine  GMTRY  or  see  Reference  (5) 
possible  combinations  of  these  para- 
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motors.  Next  read  in  the  boundary  altitudes  for  the  FASCODE 
layers  (required).  Format  (8K10.3) 


CAHD  3:  VI,  V 2 

(?F10. 3) 

VI, VP  INITIAL  AND  FINAL  WAVENUMBERS  FOR  USE  IN  CALCULATING 
THF  INDEX  01-  REFRACTION  (CM-1) 

IF  MODEL  =  7,  THE  INPUT  ATMOSPHERIC  PROFILE  13  READ  IN  AFTER 
CONTRAL  CARD  1  IN  THE  FOLLOWING  FORMAT 

( HEADER ( J ) ,  I  =  1 , P )  (2A10) 

A  PO  CHARACTER  HEADER  DESCRIBING  THE  PROFILE 
X,P,T  (3F10.3) 

(DENSITY (K) ,  K  =  1,  KMAX )  (8E10.3): 

TWO  CARDS  FOR  EACH  OF  THE  IMOD  LEVELS  GIVING  TIIF- 
ALTJTUDE  (KM),  PRESSURE  (MB),  TEMPERATURE  (K),  AND 
DENSITIES  OF  THE  MOLECULAR  SPECIES  (MOLECULES  CM-3-) 
AT  EACH  LEVEL) 

For  a  horizontal  path: 

CARD  ? :  Z,  1',  T,  RANGE  (DEN(K),  K  =  1,  KMAX) 

( 4F10 . 3 , /, ( 8E10 . 3) ) 

Z  ALTITUDE  (KM) 

P  PRESSURE  (MB) 

T  TEMPERATIJAHE  (K) 

RANGE  PATH  LENGTH  (KM) 

DEN(K)  DENSITY  OF  THE  K'TH  MOLECULAR  SPECIES 
(MOLEUCLKS  CM-3) 

For  MODEL  =  1  to  7,  only  Z  and  RANGE  are  used  and  P,  T,  and 
DEN  are  interpolated. 


CARD  3:  NOT  USED 

For  model  7,  the  input  model  atmosphere  is  read  in  after  control 
card  1  as  for  a  slant  path. 

Sample  ATM  PTH  input,  control  card  sequences  are  given  as  part  of 
the  FASCODP  users  guide  in  the  Appendix. 

Verification  of  ATMI’TH  As  A  Subroutine?  Of  Fascode :  A TM PT H  as 
supplied  by  AFGL  was  capable  of  being  run  as  an  off-line  program 
whose  output,  consisting  of  atmospheric  information,  could  be 
read  directly  by  FASCOD1 .  To  verify  that  ATMPTH  has  been  suc- 
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cessfully  interfaced  with  FASCODE  as  a  subroutine,  comparison 
runs  for  horizontal,  vertical,  and  slant  (long  and  short)  paths 
were  made  using  the  latest  version  of  FAHC0D2  with  ATMPTH  as  a 
subroutine  and  an  original  version  of  FASC0D1  using  ATMPTH  as  an 
off-line  program.  The  codes  produced  identical  output  in  all 
caseo,  indicating  that  the  atmospheric  information  was  being 
properly  transferred  to  FASC0D2  from  the  newly  added  ATMPTH  sub¬ 
routine. 
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3.  ATMOSPHERIC  AEROSOL  MODELS. 


FASCODE  performs  a  line  by  line  calculation  of  molecular 
transmittance  and  radiance  through  a  layered  atmosphere.  Out¬ 
side  the  microwave  region^,  it  does  not  include  the  effects  of 
molecular  scattering  or-  aerosol  scattering  and  absorption. 

These  extinction  mechanisms  have  now  been  added  as  part  of  the 
development  of  FAR COD? .  A  new  subroutine,  A ERROL,  gives  the 
FASC0D2  user  the  option  of  choosing  any  of  the  aerosol  models 
developed  by  AFGL  for  the  LOWTRAN  code.  The  user  is  referred 
to  Reference  8  for  a  complete  description  cf  the  available  aero¬ 
sol  models.  Molecular  scattering  is  also  accounted  for  in 
AERSOL  using  an  expression  for  the  attenuation  coefficient  taken 
directly  from  LOWTRAN. 

Description  of  the  AERSOL  Subroutine:  The  AERSOL  subroutines 
AERFRP,  FREDTA,  EXABIN,  EXTDTA,  and  AGHF.XT  have  been  taken  di¬ 
rectly  from  LOWTRAN  5.  The  only  changes  that  have  been  made  to 
the  LOWTRAN  routines  have  been  the  reduction  m  size  and  number 
of  common  blocks  wherever  LOWTRAN  variaoies  were  not  used  by 


FASCODE. 

AERSOL  is  called  initially  to  generate  the  aerosol  den¬ 
sities  at  the  31!  altitudes  of  the  model  atmospheres.  During 
this  first  call  the  single  aerosol  control  card,  which  contains 
IHAZE,  ISEASN,  IVULCN,  and  VIS,  is  read.  Subroutine  AKRPRF  is 
then  called  3^  times  (once  for  each  altitude)  to  load  the  de¬ 
sired  aerosol  density  from  the  data  stored  in  subroutine  PRFDTA , 
These  densities  are  unitlesr.  and  used  to  obtain  equivalent  sea 
level  absorber  amounts.  The  chosen  profile  is  then  stored  in 
the  array  EHM(3i0  in  the  common  blocx  /PROF/.  At  this  point 
control  returns  to  FA3C0D2,  where  ATM PTH  is  called  to  generate 
the  layer  data.  EHM^^)  is  passed  to  ATM  PTH  where  it  is  handled 
•in  the  same  way  as  the  molecular  density  profiles.  The  inte¬ 
grated  aerosol  amounts  are  stored  in  the  array  AWKAER(L)  and  are 
assigned  to  the  simple  variable  name  WKAER  in  subroutine  PATH. 


lb 
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Its 


AERSOL  is  subsequently  called  once  for  each  layer, 
function  then  is  to  load  molecular  scattering  and  aerosol  scat¬ 
tering  and  absorption  "effective  optical  depths"  into  the 
FASC0D2  array  ABSRB(J),  over  the  frequency  range  of  V1ABS  to 
V2ARS  in  increments  of  DVABS,  as  specified  by  FASC0D2.  Note 
that  this  array  was  previously  used  by  FASCODE  for  molecular 
continuum  data  such  as  H^O  and  N£.  it  now  contains  in  addition 
to  the  continuum  data,  molecular  scattering  and  aerosol  extinc¬ 
tion  which  are  also  relatively  slow  functions  of  frequency. 

The  effective  optical  depths  are  given  by  the  product  of 
the  integrated  amount  and  the  appropriate  attenuation  coeffi¬ 
cient.  For  aerosol  scattering,  or  absorption,  the  amount  WKAER 
is  known  from  ATMl'Til.  The  coefficients  are  a  function  of  the 
aerosol  model  chosen,  the  relative  humidity  (for  altitudes  be¬ 
tween  0  and  2  km),  the  altitude  region,  and  the  frequency. 

During  AERSOL* s  first  call  within  the  layer  loop  EXABIN  is  called 
to  obtain  a  set  of  altitude  and  frequency  dependent  attenuation 
coefficients  using  the  data  stored  in  subroutine  EXTDTA.  These 
are  stored  in  the  arrays  AB3C(4,4o)  and  EXTC(4,40),  correspond¬ 
ing  to  the  4  altitude  regions  and  40  wavelengths  contained  in 
VX2 ( 40 ) . 

The  frequency  loop  in  aersol  calls  AEREXT  to  obtain  the 
coefficients  at  each  frequency  by  interpolation  of  the  ABSC  and 
KXTC  arrays.  The  aerosol  optical  depths  are  then  simply: 


I.  „  -  (EXT-ABS)  *  WKAER 

A, 

ta  a  =  ABS  *  WKAER 


The  expression  for  the  molecular  scattering  optical  depth 
was  taken  directly  from  LOWTRAN  as: 

(3.1) 


m, 


=  vV(9.26Y99*  1018  -  1.07123*109*v2) 


m  - 
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For  the  case  where  the  current  FASCOD2  layer  is  not  totally 
contained  within  one  of  the  aerosol  regions,  0-2  km,  2-10  km, 

10-30  km,  or  30-100  km,  a  linear  combination  of  the  altitude  de¬ 
pendent  aerosol  attenuation  coefficients  is  used. 

The  overall  flow  of  the  aerosol  routine  is  shown  in  Figure  2. 
AEHSOL  As  A  Subroutine  Of  FASCODE:  As  mentioned  earlier,  AERSOL 
is  called  once  to  initialize  the  aerosol  density  profile  and 
then  once  for  each  layer  to  calculate  the  effective  optical 
depths  and  load  them  Into  the  ABSORB  array.  The  R1  array,  which 
eventually  has  ABSORB  merged  into  it,  now  contains  total  extinc¬ 
tion  data  rather  than  just  molecular  absorption  data.  The  cal¬ 
culation  procedure  in  FASCOD2  is  correct  for  transmittance 
but  must  be  modified  for  emission  calculations.  Emission  is 
proportional  to  absorption,  not  extinction,  therefore  the  state¬ 
ment  in  subroutine  EMIN  which  calculates  blackbody  radiation 
must  be  modified  by  the  scattering  albedo.  The  original  equation, 

BBRAD(I)  =  ( 1-TR ( I )  )  *  BB  (3-2) 

has  been  changed  to 

BBRAF(I)  =  (1-TR ( I ) )  *  (1-ALB)  *  BB.  (3-3) 

where 

k 

ALB  =  ~  . 

E 

The  user  is  referred  to  References  5  and  8  for  guidance  in 
choosing  the  proper  aerosol  model  for  his  particular  problem. 

Once  the  choice  is  made,  a  single  aerosol  control  card,  having 
the  following  format,  controls  the  aersol  routine: 


MAZE,  1SEASN,  IVULCN,  JP,  VIS 
FORMAT  (4(4x,Il),  5x,  E10.3) 


XHAZE 


=  0  no  aerosol  attenuation  included  in  the  calculation. 

=  1  RURAL  extinction,  23  km  VIS. 

=  2  RURAL  extinction,  5  km  VIS. 

=  3  MARITIME  extinction,  ?3  km  VIS. 

=  A  MARITIME  extinction,  5  km  VIE. 

=  5  URBAN  extinction,  0  km  VIS. 

=  6  TROPOSPHERTC  extinction,  50  km  VIS. 

=  7  USER  DEFINED  extinction,  23  km  VIS.  'Read  into 

the  program  immediately  after  CARD1.  Refer  to  the 

S 

main  program  LOWEM  in  Appendix  A  for  the  input 
format  of  the  coefficients). 

-  8  F0G1  (Advection  Fog)  extinction,  0.2  km  VIS. 

=  9  F0G2  (Radiation  Fog)  extinction,  0.5  km  VIS. 


ISEASN  =  0  season  determined  by  the  value  of  MODEL; 

SPRING-SUMMER  for  MODEL  =  0,1, 2, A, 6, 7 
FALL-WINTER  for  MODEL  =  3,5 
-  1  SPRING-SUMMER 
=  2  FALL-WINTER 

IVULCN  =0,1  BACKGROUND  STRATOSPHERIC  profile  and  extinction 
=  2  MODERATE  VOLCANIC  profile  and 
AGED  VOLCANIC  extinction 
=  3  HIGH  VOLCANIC  profile  and 
FRESH  VOLCANIC  extinction 
=  A  HIGH  VOLCANIC  profile  and 
AGED  VOLCANIC  extinction 
=  5  MODERATE  VOLCANIC  profile  and 
FRESH  VOLCANIC  extinction 


JP  =  print  option  parameter  to  be  used  in  the  final  version 
of  the  code;  JP  may  be  ignored  at  this  time. 

VIS  =  meteorological  range  (km) 

(when  specified,  supercedes  default  value  set  by 
THAZF.)  . 


Verification  of  The  ALR5QL  Routine:  To  verify  that  the  aerosol 
routines  function  properly  as  subroutines  of  FASC0D2,  comparisons 
were  made  with  LOWTRAN  results.  LOWTRAN  5  gives  separate  aero¬ 
sol  contributions  to  transmittance.  FASC0D2  does  not  do  this, 

so  statements  were  added  to  artificially  suppress  molecular  ab- 

20 

sorption  by  dividing  molecular  absorber  amounts  by  10  .  FASC0D2 

was  then  run  for  transmittance,  and  the  results  checked  with 
LOWTRAN  5  results  for  aerosol  absorber  amounts,  attenuation  co¬ 
efficients,  and  transmittance.  Vertical  paths  were  used  to  elim¬ 
inate  any  differences  caused  by  the  different  geometry  routines 
now  used  in  these  two  codes. 

The  comparisons  made  are  convincing  evidence  that  FASCOD2 
handles  the  aerosol  models  and  calculations  correctly.  The  im¬ 
proved  geometry  routine  now  used  in  FASC0D2  uses  a.  slightly  dif¬ 
ferent  humidity  calculation  for  the  0-2  km  region  than  used  in 
LOWTRAN  5.  This  causes  some  very  small  differences  in  the  num¬ 
bers  returned  by  the  two  codes  for  aerosols  only.  Since  the  im¬ 
proved  geometry  routine  is  being  added  to  LOWTRAN  by  AFGL,  a 
common  humidity  claculation  will  be  agreed  upon  and  put  into 
both  codes. 
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NON-LTE  TRANSMISSION  AND  RADIANCE 


Introduction  And  Formalism:  The  KASCObK  algorithm  constructs 
molecular  spectral  absorption  functions  using  an  approximate 
Voigt  line-shape  profile  for  each  active  transition  between  mo¬ 
lecular  vibrational-rotational  (VR)  states.  The  line-strength 
and  half-width  associated  with  each  transition  are  evaluated  at 
the  density-weighted  average  temperature  and  pressure  of  each  of 
the  homogeneous  layers  that  together  approximate  the  inhomogene¬ 
ous  atmosphere.  The  line  broadening  is  principally  determined 
by  collisions  in  the  higher-pressure  lower  atmosphere  (Lorentz 
region)  and  by  the  thermal  equilibrium  velocity  distribution  In 
the  upper  atmosphere  (Doppler  region).  Line-strengths,  on  the 
other  hand,  are  mainly  determined  by  the  populations  cf  the  two 
VR  states  which  each  transition  connects. 

Under  conditions  of  local  thermodynamic  equilibrium  ( LTE ) , 
collisions  cause  r-he  molecular  VR  states  to  have  a  Boitzman  pop¬ 
ulation  distribution.  In  this  case,  line-strengths  can  be  cal¬ 
culated  in  terms  of  the  strengths  at  a  reference  temperature  T  , 
since  one  knows  how  the  upper  and  lower  state  populations  vary 
with  temperature.  FASCODE  performs  this  calculation  together 
with  the  half-width  calculation  as  explained  in  Ref.  1.  Ref¬ 
erence  line-strengths  and  half-widths,  as  given  by  the  line 
parameters  compilation J  at  a  reference  temperature  (296°K)  and 
pressure  (1013mb),  are  modified  to  give  the  actual  strength  and 
both  Doppler  and  Lorentz  widths  according  to  the  average  layer 
temperature  and  pressure.  Under  LTE  conditions,  the  layer  trans 
mission  spectral  function  is  calculated  line-by-line  by  H1RAC 
and  LBLF4 .  The  layer  radiance  is  then  directly  calculated  from 
the  optical  depth  for  molecular  absorption  using  the  Planck 
function,  which  represents  thermal  equilibrium  population  ratios 
evaluated  at  the  average  temperature  of  the  layer  T. 

The  transmittance  and  radiance  in  the  lower  atmosphere  can 
be  represented  quite  adequately  in  terms  of  the  LTE  model  out- 
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lined  above.  However,  at  altitudes  above  30-50  km,  the  lower 
gas  density  enhances  the  characteristic  time  between  collisions 
and  slows  the  rate  of  VR  state  'mixing'.  The  molecular  VR 
state  populations  tend  to  be  driven  from  equilibrium  by  the  local 
radiation  field.  Although  collisions  are  not  sufficient  to  bring 
about  a  thermal  population  distribution  among  the  vibrational 
states,  both  translational  and  rotational  degrees  of  freedom 
largely  remain  in  thermal  equilibrium  at  the  kinetic  temperature 
T.  Thus,  the  line  half-widths  and  line-shapes  are  not  signifi¬ 
cantly  altered  in  the  higher  altitude  non-LTE  case. 

A  computer  code  capable  of  modeling  the  principal  non- 
equilibrium  vibrational  populations  of  CO ^ >  0^  and  NO  has 

been  previously  developed^.  We  are  adding  to  PASCODE  the  capa¬ 
bility  to  use  this  high-altitude  vibrational  state  population 
data  in  order  to  allow  for  a  more  realistic  calculation  of  both 
transmission  arid  radiance  in  this  region.  Minor  modifications 
can  be  made  to  the  basic  PASCODE  algorithm  to  allow  for  the  line- 
by-line'  construction  of  the  spectral  absorptance  function  for 
non-equilibrium  vibrational  populations.  Several  more  extensive 
changes  will  allow  the  code  to  handle  non-LTE  radiance. 

In  the  non-LTE  case,  the  general  expression  for  the  strength 
of  a  single  line  in  terms  of  the  line-strength  at  the  reference 
temperature  is 


—  Q  ( rn  ’S 


(*».!) 


where  the  p's  are  upper  and  lower  single  VR  state  populations 
for  the  homogeneous  layer  and  the  superscript  'o’  refers  to  equi¬ 
librium  conditions  at  the  reference  temperature  T  .  The  total 

optical  depth  t  is  obtained  by  summing  the  single  line  absorp- 
v  . 

tion  coefficients  K  ,  weighted  by  the  absorber  amounts  W^,  over 
‘all  of  the  active  lines  for  all  molecules: 


(A. 2) 
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The  spectral  radiance  per  line  R*  for  a  thin  layer  of  gas  (no 
self-absorption)  is  related  to  the  absorption  coefficient  for 
that  line  by 


(4.3) 


giving  a  total  emission  R^  =  X  R^W^  for  all  11 
layer.  Note  that  at  equilibrium,  the  ratio  of 
sorption  reduces  to  the  Planck  function 


nos  in  a  thin 
emission  to  ab- 


(4.4) 


for  each  line.  The  Planck  function  then  also  similarly  relates 
total  LTE  emission  and  absorption.  This  general  relationship 
does  not.  apply  in  the  non-LTE  case. 

It  will  be  assumed  that  the  line-shapes  for  absorption  and 
radiance  in  the  non-LTE  problem  are  the  same  as  they  would  be  at 
equilibrium,  so  the  K^'  and  differ  from  their  equilibrium  coun¬ 
terparts  only  by  factors  which  depend  on  VR  single  state  popula¬ 
tions.  The  actual  populations  p^  and  pu  for  a  particular  layer 
can  be  calculated  from  the  vibrational  state  population  data  of 
the  Degges  model,  assuming  the  rotational  sub-states  follow  the 

Boltzman  distribution.  Reference  single  state  populations  p° 

o  ^ 

and  p  can  also  be  calculated  for  the  thermal  equilibrium  dis- 

u  i  -i 

tribution  at  T  .  Therefore,  both  K  and  R  can  be  calculated 
o  '  v  v 

line-by-line  and  properly  summed  to  give  the  total  optical  depth 
and  radiance. 

The  model  we  then  have  is  one  containing  an  equilibrium  gas 
A  and  a  non-equilibrium  gas  B.  At  lower  altitudes,  only  gas  A 
is  active.  Above  30-50  km.,  the  majority  of  lines  are  still  in 
equilibrium,  but  some  lines  must  be  included  in  gas  B  rather 
than  A.  The  total  optical  depth  is 
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A 


+  c 


(4.5) 


B 

v  » 


where  both  and  x  ^  must  be  calculated  line-by-line, 
total  radiance  for  a  thin  layer  is 


The 


K 


v 


+  H 


B 


(4.6) 


where  F3v(T),  and  is  calculated  line-by-line.  The 

total  radiant  intensity  including  self-absorption  for  a  finite 
homogeneous  layer  is 


RAD 
J  v 


(4.7) 


A  Preliminary  Outline  Of  Planned  Code  Modifications  For  Non-LTE 
Transmission  and  Radiance:  We  describe  somewhat  schematically 
in  this  section  a  series  of  code  changes  which  will  serve  to  pro¬ 
vide  BARCODE  with  the  capability  to  use  non-equilibrium  vibra¬ 
tional  level  populations  for  improved  radiance  and  transmission 
calculations.  Every  effort  will  be  made  to  have  the  code  func¬ 
tion  as  it  does  presently  for  calculations  under  LTE  conditions. 
This  will  be  accomplished  by  constructing  new  routines  when  nec¬ 
essary  rattier  than  modifying  existing  ones. 

The  non-LTE  calculation  will  require  two  new  data  bases. 

9 

Vibrational  level  population  data  of  Degges  and  Smith  will  be 
read  and  stored  for  those  levels  which  are  not  at  thermal  equi- 
liorium.  The  data  should  cover  an  altitude  range  appropriate 
to  FASCODE  and  it  can  be  expected  to  vary  with  the  solar  angle 
under  daylight  conditions  and  with  the  change  to  night,  condi¬ 
tions.  The  second  data  base  will  consist  of  a  line  parameters 
file  which  will  have  an  expanded  format  to  include  a  vibrational 
transition  ID  or  some  additional  labels  to  facilitate  matching 
lines  with  corresponding  non-LTE  levels. 

Along  with  the  new  data  bases,  several  new  routines  are 
being  created.  A  now  routine  will  read  the  non-LTE  line  data 
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and  vibrational  level  populations.  It.  will  prepare  two  new  sets 
of  line-strengths  for  use  in  the  non-LTK  ubsorptnnce  and  radi¬ 
ance  calculations.  Another  routine  will  inter  e.a  I  m  1  at  e  tiie  to¬ 
tal  radiance  for  each  completed  layer. 

The  complete  program  will  function  an  follows.  At  lower 
altitudes,  H1RAC1  will  ignore  the  non-LTK  data  in  the  ro vised 
line  file  and  proceed  without  modification.  The  non-LTK  routines 
will  be  triggered  above  a  certain  altitude.  At  this  altitude 
and  above,  H1RAC1  will  calculate  absorptance  ] Ine-by-line  using 
the  appropriate  non-LTE  line-strengths  where  necessary.  A  sec¬ 
ond  line-by-line  calculation  will  be  performed  and  the  results 
stored  on  a  second  file  fo  '  later  use  in  the  radiance  calcula¬ 
tion.  This  second  calculation  can  be  done  by  II1RAC1,  using 
specially  modified  line-strengths.  The  radiance  routines  will 
be  modified  to  give  the  pro}  •  total  layer-  transmittance  and 
radiance  using  the  two  storeu  spectral,  functions.  The  single 
layer  results  will  then  be  merged  with  those  of  previous  layers 
using  the  current  algorithm. 
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5.  A  MODEL  POE  THE  CO,,  CONTINUUM 

Tne  line— shapes  of  CO.,  and  most  other  molecules  me  usually 
taken  to  be  Lorentzian  In  the  high  pros sure ,  co]  1  1  si  on— broadened 
limit.  The  two-parameter  Lorentzian  shape  function  is  completely 
specified  by  a  strength  E,  which  is  the  integral  of  the  shape 
function,  and  the  line  width  at  halt  maximum  a.  The  experimental 
evidence  suggests  that  the  actual  shape  function  for  C02  in  ap¬ 
proximately  Lorentzian  within  about  15  half-widths  of  the  line 

center  ,  but  becomes  sharply  sub-Lorentzi an  at  greater  distances 

1112 

from  the  central  peak  5  .  This  effect  can  be  understood  in 

light  of  the  limitations  of  eollisional  models  that  give  the 
Lorentzian  line-shape.  A  simple  impact  approximation  which  as¬ 
sumes  that  collisions  occur  instantaneously  gives  the  Lorentzian 
shape  with  a  half-width  a  'v  x  ^ ,  where  t  is  the  average  time  be- 

1  a  \  lj 

tween  collisions  *  .  The  large  frenuenev  separations  (v  -  v  ) 

~  O' 

typical  of  the  far  wings  correspond  to  very  short  time  scales 
which  may  be  shorter  than  the  actual  average  collision  duration. 
Therefore,  it  is  not  surprising  that  the  impact  approximation 
and  the  Lorentzian  shape  associated  with  it  should  not  properly 
describe  behavior  in  the  far  wings. 

The  initial  approach  that  we  will  use  to  better  model  the 
real  line-shape  in  the  far  wings  is  to  replace  the  Lorentzian 
shape  function  f^  by  the  more  general  mult ipurameter  shape  func¬ 
tion  f  in  the  far  wings: 

f  =  lLe-a[!v-v0l  -  "mln]b  i  if  I I  >  vmin  ■  <5-l) 

The  parameters  a,  b,  and  v  ,  can  be  determined  from  the  exper- 

11  min  12 

imental  data  of  Winters  and  Lurch  .  Winters  has  used  this 
shape  function  to  fit  measurements  of  the  far  wings  above  the 
2400  cm-1  band  head.  The  data  of  Burch  suggests  that  the  same 
form  may  also  be  used  for  the  3800  cm-"'’  and  7000  cm-1  bands. 

We  will  make  the  further  assumption  (consistent  with  the  limited 
data  which  is  available)  that  all  of  the  atmospheric  C0o  lines 
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can  be  modelled  with  the  name  shape  function  f,  with  only  3,  and 
a  to  be  specified  for  cacti  line. 

The  modification  of  the  far  wing  C0o  line-shape  in  FAS CODE 
is  presently  in  its  preliminary  stages.  (Since  the  FASCODE  re¬ 
port  of  January  '78,  the  problem  of  water  vapor  continuum  ab- 

2 

sorption  is  being  addressed  through  work  at  AFGL  ).  The  sub- 
i.orentzian  shape  function  will  be  used  to  construct  a  function 
representing  the  CO,,  continuum,  in  o  fashion  similar  to  the  way 
the  water  vapor  continuum  has  been  modelled.  An  AFGL  routine  is 
available  which  combines  lines  to  considerably  reduce  the  compu¬ 
tational  time  required  to  construct  a  continuum  function  line- 
15 

by-line  The  problem  of  how  one  can  properly  'match'  the  pres¬ 

sure-dependent  continuum  functions  to  an  efficient  line-by-line 
calculation  of  the  other  four  functions  which  represent  the  line 
center  remains  to  be  more  adequately  addressed. 


6.  ATMOSPHERIC  LAYERING 


The  interim  code,  FASC0D2 ,  models  the  i  rihomogeneous  at¬ 
mosphere  with  a  series  of  homogeneous  layers.  The  user  jy  re¬ 
quired  to  specify.  In  addition  to  the  path  and  model  atmosphere 
parameters,  the  altitudes  where  the  layer  boundaries  are  to  be 
placed.  The  geometry  package  ATMI’Tl!  then  'traces’  the  (gener¬ 
ally)  curved  path  through  each  layer  and  computes  an  average 
pressure  P  and  temperature  T  plus  molecular  and  aerosol  amounts 
for  each  layer.  Presently,  the  user  has  no  guidelines  to  aid  in 
choosing  layer  boundaries.  To  ensure  that  the  result  is  suf¬ 
ficiently  accurate,  a  large  number  of  layers  can  be  used,  but 
only  at  the  expense  of  computational  time,  which  is  particularly 
disadvantageous  if  a  large  range  of  wave  numbers  is  to  be  stud¬ 
ied.  Another  problem  can  arise  if  too  few  layers  are  used. 

Aside  from  having  to  deal  with  a  somewhat  uncontrollable  error, 
it  is  possible  to  violate  the  requirement  that  the  ratio  jf 
average  pressures  in  consecutive  layers  be  less  than  two  (2). 

If  this  happens,  execution  will  stop  in  mid-program. 

We  have  planned  several,  kinds  of  changes  to  overcome  these 
problems.  This  work  has  not  been  completed,  but  we  will  outline 
the  planning  that  has  taken  place  in  tne  last  year.  We  plan  to 
make  code  modifications  In  several  stages.  The  first  and  most 
obvious  change  will  be  to  construct  a  routine  that  prepares  the 
set  of  layer  boundaries  based  on  one  or  two  conditions  which  the 
user  can  specify.  A  condition  suggested  by  AFGL  is  that  the  max¬ 
imum  pressure  and  temperature  variations  within  each  layer  be 
fixed  by  the  user.  This  kind  of  condition  can  guarantee  that  the 
approximation  of  homogeneity  within  each  layer  is  a  good  one. 

If  it  proves  feasible,  the  layering  routine  may  be  expanded 
to  base  the  choice  of  layer  boundaries  on  a  wider  set  of  criter¬ 
ia.  For  instance,  the  allowable  error  in  any  single  layer  should 
depend  on  what  part  of  the  total  absorption  or  radiance  is  due 
to  that  layer.  Because  this  error  can  vary  with  the  kind  of  path 


.Involved  and  with  the  wavenumber  region  of  interest,  it  is  very 
difficult  to  decide  how  the  layer  boundaries  should  be  chosen 
in  order  to  minimize  the  total  error.  Based  on  these  kind  of 
considerations,  we  will  undertake  a  study  of  the  error  which 
results  when  running  the  code,  by  using  the  same  path  and  vary¬ 
ing  the  layering  scheme.  The  ’error*  can  be  calculated  by  com¬ 
parison  with  a  reference  run  which  has  a  sufficiently  large 
number  of  layers.  A  measure  of  the  error  in  the  transmission 
for  the  frequency  interval  from  to  is,  for  example, 


(6.1) 


Wc  have  written  a  short  routine  to  evaluate  this  error  by  com¬ 
paring  a  statistical  sample  of  the  output  of  the  two  runs.  The 
error  can  be  studied  as  a  function  of  the  number  and  position 
of  layer  boundaries,  the  type  of  path  specified,  and  the  wave- 
number  region.  The  goal  of  this  parametric  study  will  be  to 
develop  a  more  adequate  set  of  criteria  for  the  specification 
of  layer  boundaries.  The  end  product  may  be  either  a  set  of 
guidelines  for  the  user,  or  a  more  sophisticated  layering  al¬ 
gorithm  which  will  function  side-by-side  with  the  geomerry 
package  * 
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APPENDIX  -  FASCODE  USERS  CHIDE 


Introduction:  In  the  course  of  adding  the  ATMP'i'H  and  AERSOL 

routines  to  FA5COD1  several  new  control  card:;  were  introduced. 

In  addition,  some  of  the  original  PA8C0D1  control  cards  have 

been  consolidated  to  simplify  the  entire  control  deck  structure. 

The  complete  input  card  sequence  is  listed  below.  The 

routine  with  which  the  card  is  associated  is  also  given. 

Card  1  from  FASC0D1  main 

XID(I),  I  =  1,7 
(10A10) 

Card  2  from  FASCOD  1  main 

IHIRAC ,  ILBLF4,  ICNTHM,  IAERSL,  1  EM IT, 

ISCAN,  IFILTR,  TPLOT,  I ATM,  MPTS,  NPTS 
( 9 ( 4X , II ) ,  25X,  215) 

Card  3  Prom  FASC0D1  main 

VI,  V2,  TBOUND,  EMISIV 
( 10E10. 3) 

Card  4  from  AERSOL,  included  only  if  aerosol  attenuation 

is  to  be  included,  i.e.,  IAERSL  =  1 

IHAZE,  ISEASN,  IVULCN,  OP,  VIS 
(4 , (4X,I1) ,  5X,  E10.3) 

Card,  5,  6,  7  and  the  layer  boundaries  pertain  to  ATM PTH, 

and  are  included  only  if  IATM  -  0 

Card  5  MODEL,  ITYPE,  IIN,  IMOD,  KMAX,  RE 
(515,  5X,  F10.4) 

The  formats  for  control  cards  6  and  7  are  different  depending 
on  whether  the  path  is  horizontal  (ITYPE  =  1)  or  slant  (ITYPE = 

2  or  3). 

For  a  slant  path- 

Card  6  HI,  HP,  ANGLE,  RANGE,  BETA,  LEN 
(5F10.4,  15) 

Next  read  in  the  boundary  altitudes  for  the  FASC0D1  layers, 
with  format  (8F10.3). 
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Card  7  VI,  V2, 

(2F10. 3) 

For  a  horizontal  path- 

Card  6  Z.  P,  T,  RANGE  (DEN(K),  K  -  1,  KMAX) 

(4F10. 3,/,  (8E10.3) 

For  a  model  =  1  to  7,  only  Z  and  RANGE  are  used  and  P,  T  and  DEN 
are  interpolated.  For  model  =  0,  P  and  T  are  included. 

Card  7  not  used 

For  model  7,  the  input  model  atmosphere  is  read  in  after  control 

card  5  as  for  a  slant  path; 

( HEADER ( 1 )  I  =  1,2)  (2A10) 

Z,P,  T  (3F10.3) 

DENSITY  (K),  K  =  1,  KMAX  (8E10.3) 

Repeat  che  last  two  cards  for  each  of  the  IMOD  levels. 

Card  8  from  the  scanning  routine;  this  card  is  included 

only  if  I  SCAN  /■  0 

HWHM,  VI,  V2,  JEMIT,  JFN,  JVAR 
(3F10.3,  3(9X,I1) 

Card  9,  10,  11,  12  are  plotting  cards,  and  are  included 
only  if  1  PLOT  ¥  0 

Card  9  NAME,  PHONE,  EXTENSION,  ID 
(3A10) 

Card  10  VI,  V2,  XS1ZE,  DEJ.V,  NUMSBX,  NOENDX,  LFILE, 

LSKIPF,  SCALE,  IP'! 

(4F10.3,  ^15,  F10 . 3 j  15) 

Card  11  YM1N,  YMAX,  YSIZE,  DELY ,  NUMSBY,  NOENDY,  IDEC, 
JEMIT,  JPLOT,  LOGPLT 
(4F10.3,  615) 

Card  12  Repeat  card  10 

If  the  user  is  supplying  layer  data  (JATM  ¥  0)  then  FILE 
7  must  contain  the  following  information  to  be  read  by  PATH. 

NLA  YUS  0  5) 

The  next  two  cards  are  repeated  for  J  layers 

PliAR(.l),  TBAR(J),  DUMMY,  ]CNTRK(J) 

(3FJ0OI,  15) 

( AM0UN1 ( 1 ,J)  1  »  1,8) 

(8E10. 3) 


3^ 


If  aerosol  attenuation  is  to  be  included  it  must  also  be 


supplied  by  repeating  the  following  card  for  .1  layers 
AWKAER(J)  (E10.3) 

Sample  Input  Sequence:  The  fol 1 ow ing  smaple  Input  card  se- 
quency  corresponds  to  a  FASCODl  run  in  the  emission  mode,  in¬ 
cluding  aerosols  for  a  vertical  path  from  0  to  ?0  Urn,  with  a 
layer  boundary  every  1  km.  Note  that  the  scanning  routine  is 
used  (ISCAN  =  1),  while  the  plotting  routine  which  is  machine 
dependent  is  not  used  (IPLOT  =  0). 


DATA  FOR-  FASC0D2  REFORTr  GROUND  TO  20  Mir  VERTICAL  rAF.ROSOl  S  ON 
HI  =  1  E4=l  CN  =  X  AE-1  EN-1  SC=1  EI=0  PL=0 
2.095E+03  2-,  1 05F1-03  O.OOOEl-OO  O.OOOE+OO  0.000E+C0 
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-1. 


Sample  Output:  The  output  has  been  truncated  to  include  only 
output  up  to  the  second  layer,  plus  the  final  layer  and  the 
results  of  the  scanning  routine.  This  should  be  sufficient  for 
a  user  to  verify  that  his  copy  of  the  code  is  functioning 
properly . 
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